Benchmarks of the full configuration interaction, Monte Carlo shell model, and 

no-core full configuration methods 
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We report no-core solutions for properties of light nuclei with three different approaches in order 
to assess the accuracy and convergence rates of each method. Full configuration interaction (FCI), 
Monte Carlo shell model (MCSM), and no core full configuration (NCFC) approaches are solved 
separately for the ground state energy and other properties of seven light nuclei using the realistic 
JISP16 nucleon-nucleon interaction. The results are consistent among the different approaches. 
The methods differ significantly in how the required computational resources scale with increasing 
particle number for a given accuracy. 
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I. INTRODUCTION AND MOTIVATION 

Ah initio approaches to nuclear structure and reac- 
tions for p-shell nuclei have advanced significantly in 
the last few years [H-i]. At the same time, fundamen- 
tal approaches to the nucleon-nucleon [NN) and three- 
nucleon (NNN) interactions, such as meson-exchange 
theory and chiral effective field theory, have yielded ma- 
jor advances Successful realistic TV TV interactions 
from inverse scattering have also emerged 0] . These ad- 
vances in microscopic nuclear theory combine to place se- 
rious demands on available computational resources for 
achieving converged properties of p-shell nuclei. In or- 
der to access a wider range of nuclei and experimental 
observables, while retaining predictive power, we require 
additional major advances in many-body methods. 

These considerations motivate us to investigate the no- 
core Monte Carlo shell model (MCSM) which has ad- 
vantageous scaling properties for accessing larger basis 
spaces and heavier nuclei. The MCSM was first intro- 
duced in Ref. [10] and we extend it here to treat systems 
without a core. In the present work we evaluate proper- 
ties of a set of p-shell nuclei using the no-core MCSM and 
compare with exact results in the same single-particle ba- 
sis from the full configuration interaction (FCI) method 
when feasible. We also compare with representative re- 
sults from the full space ab initio no core full configu- 
ration (NCFC) fu\ method. We adopt the JISP16 NN 
interaction [9] without renormalization and without any 
NNN interactions. 

For each of the three many-body methods, all A nu- 
cleons in the nucleus are treated on the same footing. 
Experimental observables are obtained from A-nucleon 
wave functions resulting from Hamiltonian diagonaliza- 
tion in the chosen many-body basis space. To perform 
the comparisons among the methods, we focus on ground 
state properties of seven nuclei as well as the properties 



of two low-lying narrow excited states. 

For each method, we adopt the harmonic oscillator 
(HO) single-particle basis. We obtain eigensolutions of 
the nuclear intrinsic Hamiltonian expressed as a super- 
position of Slater determinants in the HO basis (FCI and 
NCFC) or the total angular momentum projected and 
parity projected deformed Slater determinants (MCSM). 
Neutron and proton orbitals are treated independently. 
The resulting calculated ground state energy is a rigorous 
upper bound on the exact result at any truncation. This 
upper bound character applies to the lowest calculated 
state of each total angular momentum and parity. 

A major distinction among the methods is the defini- 
tion of the cutoff that defines the finite many-body basis 
space in which the calculations are performed. All three 
methods should approach the exact solutions as the cut- 
offs are removed. Both the MCSM and the FCI methods 
employ a cutoff in the single-particle basis A^sheii which is 
the highest shell of the symmetric three-dimensional HO 
that is included. All many-body basis states consistent 
with that cutoff are retained (FCI) or stochastically sam- 
pled (MCSM). On the other hand, the NCFC approach 
represents an extrapolation to the infinite matrix limit of 
a sequence of calculations in many-body basis spaces de- 
fined by a many-body basis cutoff iVmax, the maximum 
number of HO quanta included in a many-body basis 
state above the minimum for that nucleus. 

A further distinction among the methods emerges from 
these different truncations — the NCFC approach may, 
in principle (though this is not used in the present bench- 
mark), guarantee the factorization of the total wave 
function into an intrinsic (translationally invariant) part 
times a pure Os HO for the center-of-mass (cm.) motion 
whereas the MCSM and FCI approaches do not guaran- 
tee this factorization. The method of analysis introduced 
for the coupled cluster method [T^l implies that MCSM 
and FCI may factorize reasonably well at an optimally 
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chosen oscillator parameter so that observables may be 
evaluated with minimal influence from spurious cm. mo- 
tion effects. All other known symmetries of the intrinsic 
Hamiltonian are retained in the many-body basis by each 
method. 

The main motivation for the no-core MCSM approach 
is its superior scaling properties with increasing nucleon 
number. We estimate that, for a fixed iVshcii value and a 
fixed level of accuracy the MCSM scales as x N^~'^ 
where N/, is the number of Monte Carlo basis states gen- 
erated in the sampling and A'sp is the number of HO 
single-particle states included by A^shoii- To obtain a fixed 
accuracy with increasing nucleon number A, Nb will have 
to increase as some low power of A, estimated at 1.5 ^ 2.5 
from the results we present below. Assuming domi- 
nates the A dependence for fixed accuracy, which seems 
reasonable, we estimate that MCSM scales as A^"^^. On 
the other hand, the NCFC scales as ^i2--i4 f^^. ^ g^^gj 
Nma.x value and the maximum A^max value roughly fixes 
the accuracy of the final NCFC result. Since the MCSM 
scaling for fixed accuracy is far less dependent on the 
number of nucleons A, it will be the superior approach 
once A increases to the point where the NCFC fails to 
generate a sufficiently converged result. Nevertheless, the 
truncated calculations within NCFC will continue to pro- 
duce a valid upper bound to the exact answer. 

Since the MCSM approximates the FCI calculation by 
stochastically sampling the FCI many-body basis space, 
we provide comparisons between these two methods in 
smaller basis spaces and for lighter systems where we can 
still perform the FCI calculations. For these test prob- 
lems, we find that the MCSM provides an accurate ap- 
proximation to the FCI results. The sequence of MCSM 
results with increasing A^sheii and for heavier nuclei may 
also be compared with the sequence of results as a func- 
tion of A^max that underlie the NCFC result in order to as- 
sess convergence rates and uncertainties in extrapolated 
results. 

The outline of this paper is as follows. After the Intro- 
duction and Motivation of Sec. HI many-body basis space 
truncations and quantum many-body methods adopted 
for the benchmark in this paper are briefly described in 
Sec. nil The selections of the A^A^ interaction and nu- 
clear states are summarized in Sec. IIIII The benchmark 
comparisons are presented and discussed in Sec. |lVl The 
summary and outlook can be found in Sec.|Vl In the Ap- 
pendix we present additional details for the energy vari- 
ance and the extrapolation of the no-core MCSM results 
to the FCI basis. 



II. QUANTUM MANY-BODY METHODS 
ADOPTED 

A long-standing goal of nuclear physics is to obtain the 
exact solutions of a realistic Hamiltonian (i.e., one that 
describes well the few-body data) for finite nuclei and to 
compare those results with experiment where available. 



Once validated, the same methods with the same Hamil- 
tonian will be very useful for predicting properties of nu- 
clei that cannot be studied experimentally but may be 
of great importance in understanding astrophysical phe- 
nomena or for practical applications such as energy gen- 
eration. This is the physics program we aim to empower 
by developing and testing new many-body methods. 

We begin by introducing the elements that the three 
methods we study here have in common. The transla- 
tionally invariant nonrelativistic nuclear plus Coulomb 
interaction Hamiltonian is taken to consist of 

H = Trcl + VnN + VnNN -I- . . . + Vcoulomb, (1) 

where Tid is the internal ("relative") kinetic energy of 
the nucleons and the A^A^ and NNN interactions are in- 
cluded along with the Coulomb interaction between the 
protons. The Hamiltonian may include additional terms 
such as multinucleon interactions among more than three 
nucleons simultaneously and higher-order electromag- 
netic interactions such as magnetic dipole-dipole terms. 

The JISP16 A^A^ interaction adopted here produces a 
high-quality description of the A^A^ scattering data and 
the deuteron as well as a good description of a range 
of properties of light nuclei 'iTl. For the present effort 
we neglect all other interaction terms such as the NNN, 
higher-body strong interactions and the Coulomb inter- 
action though the three methods are capable of includ- 
ing them. These additional terms will be required for 
precision descriptions of nuclear properties but are not 
expected to alter the conclusions from our benchmarks 
here. 

All calculations are performed in an M-scheme basis 
where the many-body basis states are constructed with 
good total magnetic projection M. The MCSM projects 
out states of fixed total angular momentum and parity 
J'^. The basis states used in the FCI and NCFC calcula- 
tions are constructed with a fixed parity (as well as fixed 
M). The eigensolutions of the FCI and NCFC methods 
will also possess good J up to numerical errors. Eval- 
uating the value of J for any eigensolution serves as a 
crosscheck on the precision of the calculations. 

In all applications here, we seek to obtain only the low- 
est few eigenvalues and eigenfunctions. For the NCFC 
and the FCI calculations we employ the code "Many- 
Fermion Dynamics - nuclear" or "MFDn" [l^ which 
has been optimized for leadership-class parallel comput- 
ers For the MCSM calculations, we employ a new 
MCSM code that runs efficiently on parallel comput- 
ers [151. 

All solutions will have a dependence on the cutoff (ei- 
ther Ashou for FCI and MCSM or A^max for truncated 
NCFC) and dependence on the HO energy hui. The 
MCSM results also depend, in principle, on the num- 
ber of Monte Carlo basis states A";, and we employ an 
extrapolation based on energy-variance to estimate the 
Aft-independent solution. The degree to which we obtain 
results independent of the cutoff and of the HO energy 
is a measure of the convergence of the results — fully 
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converged results are independent of all basis space pa- 
rameters. 



A. Many-body basis space truncations 

The methods we investigate employ one of two differ- 
ent truncation schemes as mentioned above. The MCSM 
and FCI employ an A'shcU cutoff while the NCFC em- 
ploys TVinax to define the finite basis spaces in which 
the Hamiltonian is evaluated and diagonalized. We work 
in a neutron-proton scheme rather than a basis of good 
isospin. We now discuss some additional features of those 
truncation schemes. 



1 ■ A'shcU 

For the MCSM and FCI methods, all single-particle 
states for neutrons and protons in HO shells up to and 
including A'shoii are included (A'sheii = 1 for the low- 
est shell). Then, all many-body states consistent with 
that cutoff and the selected symmetries are enumerated. 
Thus, for example, we include basis states where all nu- 
cleons occupy the highest HO shell if that shell can ac- 
commodate all of them. Table |T] presents many-body ba- 
sis space dimensions in the M scheme and J scheme over 
a range of iVsheii values for the nuclei we investigate. We 
also include ^^O for illustrative purposes. An FCI cal- 
culation involves evaluating the Hamiltonian with that 
dimension and diagonalizing it — at least to obtain the 
low-lying solutions of interest. 

ax 

For the NCFC method, we employ the many-body 
A'niax truncation where we enumerate all many-body 
states, with the selected symmetries, possessing total HO 
quanta less than or equal Nj^ax above the lowest allowed 
configuration for that nucleus. Each single-particle state 
in a basis state contributes 2n + l to the total HO quanta 
(n is the radial quantum number and I is the orbital an- 
gular momentum quantum number) for that basis state 
and then the minimum sum for that nucleus is subtracted 
to give the total quanta above the minimum for that ba- 
sis state. The basis space for each nucleus begins with 
-^max = and increases in units of 2 for the natural par- 
ity states. Odd values of A^max cover the unnatural parity 
states. Thus, for example, we include basis states where 
one nucleon occupies the highest HO shell accessed. Ta- 
ble |ll] presents many-body basis space dimensions in the 
M scheme over a range of Amax values for the nuclei 
we investigate, again with ^^O added for illustrative pur- 
poses. A no-core shell model (NCSM) calculation in- 
volves evaluating the Hamiltonian with that dimension 
and diagonalizing it - at least to obtain the low-lying 
solutions of interest. A sequence with increasing A'max 



of NCSM calculations using a Hamiltonian defined for 
the infinite basis will converge from above to the exact 
solution. The NCFC approach uses that sequence to ex- 
trapolate to the infinite basis limit. 

B. FCI 

An FCI calculation involves solving the Hamiltonian 
eigenvalue problem in a many-body basis space with the 
A"sheii truncation described above. We have performed 
sets of these calculations in the present effort to pro- 
vide the exact results for comparison with the MCSM 
approach and to compare with the A^max-truncated re- 
sults of the NCFC approach. For the FCI results reported 
here, we employ the M-scheme basis whose dimensions 
are indicated in Table |T] and use the Lanczos algorithm 
in a manner similar to a NCSM calculation. Unlike the 
NCFC approach, we do not perform an extrapolation to 
the infinite matrix limit of the FCI results as a function 

of AshoU- 

C. MCSM 

The MCSM approach [13, [H proceeds through a 
sequence of diagonalization steps within the Hilbert 
subspace spanned by the selected importance-truncated 
bases, beginning with, in principle, any initial trial solu- 
tion for the system. Until now, the deformed Hartree- 
Fock (Hartree-Fock-Bogoliubov) states in the HO single- 
particle basis defined by the A'sheii cutoff have been 
adopted as an initial state for the shell-model calculations 
with a core in light (medium-heavy) nuclei. These de- 
formed single-particle states in a canonical basis are con- 
structed as a linear combination of spherical HO single- 
particle states up to and including those in the A'sheii cut- 
off. One then stochastically samples all possible many- 
body basis states around the mean field solutions with 
the aid of the auxiliary fields and diagonalizes the Hamil- 
tonian matrix within the subspace spanned by these 
bases. An accept/reject process of a stochastically sam- 
pled basis is performed by minimizing the energy vari- 
ationally, not by the importance sampling in quantum 
Monte Carlo methods implemented by the Metropolis al- 
gorithm. The MCSM is thus not the usual "quantum" 
Monte Carlo, but can evade the so-called negative sign 
problem, which is the fundamental issue that cannot be 
avoided in quantum Monte Carlo methods. 

In the MCSM, a many-body state j^'-^"*^) is con- 
structed from the linear combination of non-orthogonal 
angular-momentum (J) and parity (tt) projected de- 
formed Slater determinants |<I>) with good total angular 
momentum projection (M) as a stochastically selected 
basis, 

l*"''') = E (2) 
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7 


M scheme 


*He 


98 


3.06 


X 


lO'' 


3.98 X 10" 


3.14 X lO'' 


1.77 X 10" 


7.84 X 10" 


«He 


216 


6.51 


X 


10* 


3.86 X 10"^ 


9.80 X 10^ 


1.45 X 10^ 


1.47 X 10" 


"Li 


293 


8.59 


X 


lo-* 


5.08 X 10® 


1.29 X 10** 


1.91 X 10^ 


1.94 X 10" 


'Li 


400 


3.60 


X 


lO'^ 


4.51 X 10^ 


2.05 X lO'' 


4.91 X 10" 


7.50 X 10" 


^Be 


518 


1.47 
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10"^ 


3.96 X 10^ 


3.24 X 10" 


1.26 X 10^^ 


2.91 X 10" 
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1.34 
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10^ 


1.82 X 10^" 


5.02 X 10" 


5.22 X 10" 


2.78 X 10" 
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8.22 
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10^ 


5.87 X 10" 


5.50 X 10" 


1.54 X 10" 


1.90 X 10" 


16q 
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8.12 
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10® 


2.10 X 10" 


2.51 X 10^* 


5.32 X 10^^ 


3.59 X 10^* 


J scheme 


*He 


20 


2.72 


X 


10^ 


2.10 X 10-* 


1.12 X lO'' 


4.58 X 10* 


1.54 X 10'' 


«He 


35 


3.93 
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10^ 


1.37 X 10^ 


2.35 X 10® 


2.52 X 10^ 


1.93 X 10** 
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1.42 
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10* 


5.19 X 10^ 


9.05 X 10® 


9.79 X lO'' 


7.57 X 10** 




89 


3.63 
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10* 


2.73 X 10® 


8.40 X 10^ 


1.46 X 10^ 


1.69 X 10" 


^Be 


70 


6.89 
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10* 


1.08 X 10^ 


5.92 X 10** 


1.66 X 10" 


2.90 X 10" 




43 


2.20 


X 


10^ 


1.21 X 10^ 


2.21 X 10" 


1.65 X 10" 


6.67 X 10" 


12^ 


20 


2.94 


X 


10® 


1.14 X 10^" 


6.94 X 10^2 


1.38 X 10" 


1.28 X 10" 


16o 


1 


2.54 


X 


10^ 


3.26 X 10^^ 


2.46 X 10" 


3.66 X 10" 


1.84 X 10^^ 



TABLE I: Dimensions of the M scheme (top) and J scheme (bottom) many-body basis spaces for selected nuclei with the Afshoii 
truncation. The dimensions are for the natural parity states with M and J taken to be the lowest allowed value (Af — for 
even nuclei except for ®Li and ^°B where M — 1; M = 1/2 for odd nuclei, and similarly for J). 
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lO'^ 


6.15 
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1.24 X 10^ 
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5.10 
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1.44 
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10*^ 


2.22 
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10® 


2.33 
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10^ 
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1.22 X 10^ 
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1.16 
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1.60 
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1.77 
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3.26 
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4.06 X 10" 
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TABLE IL Dimensions of the A/-scheme many-body basis spaces for selected nuclei with A^'max truncation. The dimensions are 
for the natural parity states with M taken to be the lowest allowed value {M = for even nuclei except for ®Li and ^®B where 
M = 1; M = 1/2 for odd nuclei). The sequence of dimensions for unnatural parity states (odd values of Amax) lie intermediate 
the neighboring natural parity dimensions. 



where the angular-momentum and parity projected basis, 
\<^'^^')= 9KPi,KP^\ct^), (3) 

K=-.} 

and the deformed Slater determinant, 

A 

i0)=n«ii-)' (4) 

i=i 

with the vacuum |— ) and the creation operator a\ = 

'n!a=i'^^ai- The coefficient Dai is stochastically sam- 
pled by the auxihary-field Monte Carlo technique around 
the Hartree-Fock solutions. 

With increasing Monte Carlo basis dimension, the 
ground state energy of a MCSM calculation converges 
from above to the exact value — the value that would 
be obtained by diagonalization of the corresponding FCI 
basis space. The energy, therefore, always gives the vari- 
ational upper bound in this framework. 



An exploratory no-core MCSM investigation of the 
proof-of-the principle type has been done for the low- 
lying states of the Berylium isotopes by applying the ex- 
istingMCSM algorithm with a core to a no-core prob- 
lem [22l. Recent improvements on the MCSM algorithm 
have enabled significantly larger calculations [3) fl6| . 
In addition energy variance extrapolation methods have 
been introduced and tested in order to obtain precise 
results at each A^sheii cutoff We adopt these im- 

provements in the present work, and extend our earlier 
investigations Ts']. A similar work by the hybrid multi- 
determinant method is also proposed ,19i] . 

D. NCFC 

The NCFC approach flll| aims to achieve the solution 
of the nuclear many-body problem by diagonalization in a 
sufficiently large basis space that converged energies are 
accessed — either directly or by simple extrapolation. 
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snapshot comparison FCI MCSM NCFC 

cm. motion approx. approx. exact 

Spectra OK some OK 

wfns — ^ observables / / / 

Matrix dimension < 10^" < 10^° < 10^° 

Scaling with A^^-^o ^3^5 ^i2~i4 

No. parallel cores 10^ 10^ 10^ 

Comp'l bottleneck Memory CPU time Memory 



TABLE III: Overview of the current features of the three 
no-core many-body methods employed in this work. The es- 
timates of the scaling with the number of nucleons A are very 
crude and based on applications to light (p-shell) nuclei for a 
fixed accuracy. The last two lines in the table present overall 
characteristic features of the codes used for this work. 

Convergence is assessed in the two-dimensional parame- 
ter space of the basis space {hoj, A'max) and is defined as 
independence of both parameters within estimated un- 
certainties. Each observable is studied independently to 
obtain its converged value and its assessed uncertainty. 

The NCFC is both related to and distinct from the 
NCSM [5] that features a finite matrix truncation and an 
effective Hamiltonian renormalized to that finite space. 
The A^max regulator appears in both the NCFC, where 
it is taken to infinity, and in the NCSM, where it also 
appears in the definition of the effective Hamiltonian. In 
both the NCFC and the NCSM, the iVmax cutoff in the 
HO basis is needed to preserve Galilean invariance — to 
factorize all solutions into a product of internal motion 
and cm. motion components. With A^max as the regu- 
lator, both the NCFC and the NCSM are distinguished 
from the FCI and the MCSM approaches where such fac- 
torization is not guaranteed but may be approximately 
valid . In the NCFC approach the ground state ener- 
gies at any finite truncation are strict upper bounds, and 
converge monotonically to the exact result. This facili- 
tates a simple extrapolation. The approach to the exact 
result is in general not monotonic in the NCSM. 



E. Snapshot comparison 

For convenience, we present simple comparisons among 
the methods we employ in Fig. [1] and in Table IIIII In 
Fig. [T]we show, for the specific example of ^^C, how the 
many-body basis spaces both overlap and differ from each 
other as a function of increasing cutoff. To indicate an 
area of complete overlap, the red curve in Fig. [T] borders 
the A^shoU space included as a function of increasing A^max- 
On the other hand, the blue curve borders the region of 
A^max space included as a function of increasing A^shoii- 

Table IIIII presents a simple set of comparisons and 
contrasts between the methods. Since A'sheii and A'max 
roughly signify respective accuracies of the methods we 
hold them fixed to facilitate these comparisons. We em- 
phasize that these are the current features and limita- 
tions of these approaches. Additional developments un- 




2 4 6 8 10 12 14 16 18 
N^3, [NCFC/NCSM] 



FIG. 1: (Color online) Overview of the basis spaces covered 
with the many-body methods discussed here for the case of 
^^C. A'niax is defined as the number of oscillator quanta above 
lowest possible number of quanta. A'aheii is the number of 
oscillator shells counting the Os shell as the first shell. The 
MCSM incorporates an FCI space. That is, all single-particle 
states in the included shells are available to all particles with- 
out additional restrictions except for symmetry constraints. 



derway are aimed at improving each method, especially 
the MCSM and NCFC approaches. 

The main advantage of the MCSM approach is that, 
at fixed A'shcii, the increase in computational needs with 
increasing nucleon number ("scaling" in Table [Hi]) of the 
MCSM approach is much slower than that of the FCI ap- 
poach. In addition, the increase in computational needs 
with A of the MCSM approach at fixed A'shcii is signif- 
icantly slower than that of NCFC at fixed A^max- Note 
that the MCSM algorithm is CPU bound, and may be 
suitable for implementation on general purpose graphical 
processing units (GPGPUs). 

Before progressing to the detailed comparisons among 
the results from the methods we investigate here, it 
is worth noting that there are additional efforts aimed 
at accelerating the convergence of ab initio no-core 
many-body methods using basis function techniques. 
The "Importance- Truncated" no-core shell model (IT- 
NCSM) ^] attempts to estimate the contributions of 
the many-body configurations above the A'max cutoff us- 
ing sequences of perturbative contributions to the energy 
of low- lying states. The symmetry-adapted no-core shell 
model (SA-NCSM) [2lj aims to augment the basis space 
above the A^max cutoff by adding basis states of selected 
symmetry character that are preferred by low-lying nu- 
clear collective motion. Both methods are producing im- 
pressive results. It remains to be seen which method, 
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among the many under investigation, will be more effi- 
cient and for which systems and which observables. Out- 
standing challenges include the fully microscopic descrip- 
tion of clustering phenomena and extensions to ab initio 
nuclear reaction theory. 



III. SELECTIONS OF INGREDIENTS 

We have outlined above the many-body methods se- 
lected for these benchmark comparisons (FCI, MCSM, 
and NCFC). All results are obtained in a HO basis of 
single-particle states characterized by the oscillator en- 
ergy hu in MeV and the cutoff of the basis space (A'shcii 
or iVmax) defined above. We adopt the JISP16 NN inter- 
action [9] without renormalizing it to a lower momentum 
scale and wc neglect Coulomb and all other interactions. 
The contributions of spurious cm. excitation are not 
discussed here in any detail. Such contributions are ab- 
sent in conventional NCFC results for ground state ob- 
servables where we would include the standard Lagrange 
multiplier term that constrains the cm. motion to the 
Os HO state. However, for the present benchmark com- 
parisons, we have dropped the Lagrange multiplier term 
in the Hamiltonian for simplicity. We do not expect that 
spurious cm. effects play a significant role in our bench- 
mark comparisons. 



A. Interaction 

The JISP16 NN interaction is determined by inverse 
scattering techniques from the np phase shifts and is 
taken to be charge independent. JISP16 is available in a 
relative HO basis and can be written as a sum over 
partial waves 

Vnn = T's.j.T J2 \^^^) ^nfJe' {r^'^'l (5) 

S,J,T n,i,n',l' 

where tiuj = AO MeV and J ^£ + s. The HO basis state 
of relative motion is signified by \n£) and the projector 
onto the specified channel is represented by Vs,j,t- A 
small number of coefficients {Af^'f^'^i>} are sufficient to 
describe the phase shifts in each partial wave. Note that 
the JISP16 interaction is nonlocal and its off-shell prop- 
erties have been tuned by phase-shift equivalent transfor- 
mations to produce good properties of light nuclei. For 
example, JISP16 is tuned in the ^Si —^Di channel to give 
a high precision description of the deuteron's properties. 
Other channels are tuned to provide good descriptions of 
■^H binding, the low-lying spectra of ^Li and the bind- 
ing energy of ^^O. With these ofi-shell tunings to nuclei 
with A > 3 one may view JISP16 as simulating, to some 
approximation, what would appear as NNN interaction 
contributions (as well as higher-body interactions) in al- 
ternative formulations of the nuclear Hamiltonian. 



B. Nuclear states evaluated 

For this benchmark process, we select nine states of 
light nuclei that includes seven ground states and two 
excited states; ''He (0+), ^He (0+), ^Li (1+), ^Li (1/2", 
3/2"), *Be (0+), i°B (1+, 3+), and (0+). We com- 
pare results for the energy, the point-particle root-mean- 
square (rms) matter radius, and the electric quadrupole 
and magnetic dipole moments. 

Our goal here is to compare the methods at fixed fi- 
nite cutoffs. To achieve convergence of the quantities we 
evaluate will require a much larger effort than the present 
undertaking. For the benchmark process, we simply pro- 
ceed through a sequence of cutoffs for each state and each 
method and obtain results as a function of the oscillator 
energy, huj. Then, since all our methods retain the varia- 
tional principle, we select the optimal huj that minimizes 
the energy for that state and basis space cutoff. We com- 
pare the observables for that optimal hu. 

The MCSM results are compared with those of FCI 
which gives the exact results in the chosen single-particle 
model space. The FCI results are obtained by the MFDn 
code [ij, Uland the MCSM results by the newly devel- 
oped code [isl. [l6| . Note that the FCI results are not 
available for all the cases presented here due to compu- 
tational limitations of the FCI approach as indicated by 
the "Matrix dimension" entry in Table IIIII 



IV. BENCHMARK COMPARISONS 
A. Results for energies 

We present the energies obtained by MCSM and FCI 
in Fig. [2] through a sequence of A^shoU truncations. For 
the A = 4 and 6 systems we obtain results with both 
MCSM and FCI through Nshdi = 5. For the systems with 
A = 7 and A = 8 we obtain results with both MCSM and 
FCI through iVshcii = 4. Finally, for (1+, 3+) and 
^^C (0+) we are not able to obtain the FCI results be- 
yond A^ghcii = 3 due to computational limitations so the 
MCSM results at A"shcii = 4 represent predictions. That 
is, the FCI M-scheme matrix dimensions are 18.2 billion 
and 587 billion as shown in Table U for ^"^B and ^^C, re- 
spectively, at A^shoii — 4 and these dimensions exceed our 
current FCI capabilities. 

The MCSM and FCI nearly coincide in all cases where 
both are available. In fact, for A"shcii = 2 they are identi- 
cal to within machine precision. This can easily be under- 
stood because the dimension of the complete FCI basis in 
the J scheme is below 100 (except for the J'^ = | basis 
space of ^Li, see Table H]), the number of Monte Carlo 
states used in most of the calculations presented here.^ 



The Monte Carlo basis states do not form an orthogonal basis, 
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FIG. 2: (Color online) Comparisons of the energies between the MCSM and FCI along with the fully converged NCFC results 
where available. The NCFC result for the ^°B(1"'") state has a large uncertainty indicated by the grey band. The MCSM (FCI) 
results are shown as the solid (dashed) lines that nearly coincide where both are available. The extrapolated MCSM results 
are illustrated by bands. From top to bottom, the truncation of the model space is Afshcii = 2 (red), 3 (green), 4 (blue), and 5 
(purple). Note that the MCSM results are extrapolated by the energy variance with the second-order polynomials 17j. Also 
note that all of the results of ^"B and ^^C at A^shcii = 4 were obtained only with MCSM. 



Indeed, the MCSM results tend to become independent 
of the number of Monte Carlo states Nf, once Ni, is of the 
order of 10% to 20% of the dimension D of the underly- 
ing FCI basis in the J scheme. Also for iVsheii = 3 the 
MCSM and FCI results are virtually indistinguishable 
in Fig. [21 After the extrapolation of the MCSM results 
using the energy variance [l7[ , as discussed below and in- 
dicated by shaded regions in Fig. ^ we obtain also very 
good agreement with the FCI results (where available) 
for A^shcii > 4. 

In Fig. [5] we also present the NCFC results for the en- 
ergies of the A — A through 10 nuclei (black solid lines) 
for comparison. The NCFC results are obtained from 
calculations up through A^max = 14 for A = 4 and 6, 
and up through Ny^^x = 12 for A = 7 and 8, using an 
exponential extrapolation to the infinite basis space. In 
these cases, the extrapolation uncertainties in the fully 
converged NCFC results are less than the width of the 
black line. For A = 10 we employ results up through 



and can be overcomplete. 



-^max = 10 to obtain the NCFC results. The extrapola- 
tion of the 1+ state in to obtain the quoted NCFC 
result has a significantly greater uncertainty due to the 
occurrence of two close- lying 1+ states in the calculated 
spectrum 

In order to stimulate future comparisons with other 
many body methods, we present detailed results in ta- 
bles for selected values of A^shoU- For the energies we 
present results according to the method and the basis 
space cutoff in Tabic HVl All results are presented for the 
value of hw where that state is a minimum in that A'shcn 
basis, except for the NCFC results, which are, within 
the estimated numerical uncertainty, independent of any 
basis parameters. Here, we observe that the differences 
between the MCSM and FCI results is at most a few hun- 
dred keV for A'shoii = 3, which is why they are barely dis- 
tinguishable at the energy scale of Fig. [5] For A^shoU > 3, 
this difference can be of the order of an MeV or more. 
However, extrapolated MCSM results agree with the FCI 
results to within the estimated extrapolation error, with 
only one case in which the difference is larger than the 
estimated extrapolation error in Table IIVI That case is 
®Be at A'sheii = 4 where the uncertainty is 1 keV and the 
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Method 


Nsiiall — 2 




D 


A'shell — 3 






iVshoii = 4 






A''shoii = 5 






NCFC 


"He (0+) 


MCSM 

extrp 

FCI 


-25.956 
-25.956 


30 


20 


-27.914 
-27.914 


30 


100 


-28.737 

-28.738(1) 

-28.738 


30 


100 


-29.011 

-29.037(1) 

-29.036 


25 


50 


-29.164(2) 


•^He (0+) 


MCSM 

extrp 

FCI 


-13.343 
-13.343 


20 


35 


-19.186 

-19.196(1) 

-19.196 


20 


100 


-23.480 

-23.687(4) 

-23.684 


25 


100 


-25.080 

-26.086(76) 

-26.079 


25 


50 


-29.51(5) 


«Li (1+) 


MCSM 

extrp 

FCI 


-14.218 
-14.218 


20 


97 


-21.549 

-21.581(1) 

-21.581 


20 


100 


-26.757 

-27.166(16) 

-27.168 


25 


100 


-28.410 

-29.873(83) 

-29.893 


25 


50 


-33.22(4) 


(1/2-; 


) MCSM 
extrp 
FCI 


-14.459 
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20 
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20 
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FCI 


-29.755 
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-41.965 

-42.357(46) 

-42.338 


25 
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-52.239 
-54.89(16) 
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MCSM 
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-34.221 
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-46.263 

-4D.Dio(2/ ) 
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TABLE IV: Energies in MeV calculated for the seven ground states and two excited states within the MCSM and FCI methods 
using the JISP16 A''A'^ interaction. The entries of the MCSM indicate the MCSM results before the energy variance extrapolation, 
while the those of the "extrp" line denote the MCSM results after the extrapolations. The number of Monte Carlo vectors 
evaluated in the MCSM approach is indicated by A^;,. The value cited for fio; (units are MeV) represents the value at which 
the energy for that state reaches its minimum value in that Afshdi basis. Uncertainties in extrapolated results are quoted in 
parenthesis. 



difference is 2 keV. 

The energies converge uniformly from above as ex- 
pected with increasing iVshoU- We obtain significant in- 
creases in binding with each increment in A'shoii and this 
encourages us to develop the MCSM further in order to 
access larger A'shcii bases. At the present time, our lim- 
ited results do not indicate a pattern that we can extrap- 
olate to the infinite A'sheii limit. However, the expected 
outcomes of such extrapolations should be the NCFC re- 
sults shown in Fig. [5] and in the last column of Table HVl 
Larger iVshcU results and extrapolations to the infinite 
iVshcU limit constitute goals for future efforts since our 
main goal here is to benchmark the MCSM approach 
through the range of A^shoii values accessible by FCI and 
to compare with the fully converged NCFC where avail- 
able. 

The detailed convergence pattern for ground state en- 
ergy for "^He is shown for the FCI and NCFC methods in 
the left panel of Fig. [3] as a function of fiw and the ba- 
sis space cutoff (iVghoii for FCI; iVmax for the approach 
to NCFC). We define convergence as independence of 
both hw and the basis space cutoff. We note that FCI 
at A'shoii = 8 and the NCFC truncated at iVmax = 10 



both yield almost the same ground state energy of —29.15 
MeV, even though the dimensions are quite different: the 
full iVshcU = 8 basis space dimension of ''He is 29 031 044, 
whereas the A'max — 10 basis space dimension is only 
196 438, more than two orders of magnitude smaller. The 
NCFC result for ''He is -29.164(2) MeV. For comparison 
the MCSM results at A/'sheii = 5 (the largest MCSM space 
reported here) and hut = 25 MeV, once extrapolated with 
the energy-variance method, produces —29.037(1) MeV 
which agrees to within 1 keV of the FCI result for that 
space 

We also show the detailed convergence patterns for the 
ground state energy of ^Li in the right panel of Fig. [3] for 
FCI and NCFC at various truncations as a function of haj 
and the basis space cutoff. We note that the convergence 
trends of the A^shcii = 2 — > 5 results for ^Li shown in the 
right panel of Fig. [3] has both similarities and differences 
from the pattern for ^He seen in the left panel. Both ex- 
hibit the "U" -shaped patterns for each truncation with 
the bottom of the "U" becoming flatter as the cutoff in- 
creases. However, for ^He, the minimum with respect to 
huj for iVshcii > 1 remains at nearly a constant hu value as 
the cutoff is removed while for ^Li that minimum shifts 
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to higher values of hu. 

For ^Li, the ground state energy increments by about 
the same amount from A'shcii = 3 ^> 4 as for A'shoii = 2 — > 
3. However, there is a significant decrease in the energy 
increment for the step iVshoii = 4 — )■ 5. Furthermore, we 
observe that the energy increment for A'shcU = 4^5 
is approximately the difference between the A'shcii — 5 
result and the fully converged result given by NCFC. It 
will be valuable to extend the A'shcii cutoff further in a 
future effort in order to determine the full convergence 
pattern for ^Li. Note that the results obtained in the 
FCI Nshe\i = 5 basis space, with a dimension of 129 mil- 
lion, are very close to the iVmax — 6 results, with a basis 
space dimension of less than 0.2 million. On the other 
hand, with only 50 Monte Carlo basis states, the MCSM 
produces a ground state energy that is within 0.5 MeV 
of the FCI result at A'shcii = 5. With extrapolation, the 
ground state energy is within 20 keV of the FCI result 
with an extrapolation uncertainty of 83 keV. 

B. Convergence of MCSM calculations 

In top panel of Fig. H] we show the convergence of the 
ground state energy of ^He as function of the number 
of the Monte Carlo basis states, Nf,. At every value of 
Nf, the MCSM gives a variational upper bound for the 
energy, and as Ni, increases, the energy approaches the 
exact FCI result from above. In the bottom panel we 
show the relative difference between the MCSM and the 
FCI result, not only for the energy but also for the rms 
radius. 

Both the ground state energy and the radius of ^He 
converge to within 1% of the exact FCI result with less 
than ten Monte Carlo basis states, at least up to A^shoii — 
5. However, in general, as the number of shells increases, 
so does the number of Monte Carlo basis states that is 
needed in order to achieve a fixed level of accuracy: at 



^shoii = 4 we need about 20 basis states in order to reach 
an agreement of 0.1%, but at A^sheii = 5 we need about 
50 basis states to reach the same level of accuracy, as can 
be seen from the bottom panel of Fig. [4l 

The number of the Monte Carlo basis states needed for 
a given level of accuracy depends not only on iVshcii, but 
also on the number of nucleons A, the quantum numbers 
of the state under consideration, and the observable, as 
can be seen from Fig. [SJ 

In general, the convergence with Nj, to the exact FCI 
result starts out very rapidly, but slows down as Nt in- 
creases. The energy always converges monotonically (at 
least for the lowest states of a given spin and parity) , be- 
cause of the variational principle, but other observables 
such as the rms radii and the magnetic moments do not 
converge monotonically. On average, however, the dif- 
ference between the MCSM results and the FCI results 
decreases with increasing Nt, as one would expect. 

Furthermore, the average convergence rate with in- 
creasing Nb for different observables of a particular state 
at fixed A^shoU tends to be the same. That is, if the MCSM 
energy converges rapidly to the FCI energy, then so do 
the rms radius and magnetic moments of that state, but 
if the MCSM energy converges slowly, then the other ob- 
servables converge slowly as well, as one can see from 
Figs, m and [51 This suggests that the wave function ob- 
tained with the MCSM converges to the FCI wave func- 
tion in a systematic manner that can be measured by 
different observables. 

In Fig. [S] we show the number of Monte Carlo basis 
states Ni, that are needed in order to achieve a specified 
level of accuracy for the energy as function of A for the 
four nuclei under consideration that have a 0"*" ground 
state. Clearly, the convergence for '''He is much faster 
than for any of the other nuclei, and its convergence rate 
is (unfortunately) not a good indicator for the conver- 
gence rate that can be expected for heavier nuclei. For 
the other three nuclei we see that at iVshcii = 3 a dou- 
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FIG. 4: (Color online) The convergence of the MCSM ground 
state energy of *He to the FCI result for several A'aheii val- 
ues as function of the number of Monte Carlo basis states, 
Ni,. Top: ground state energy; Bottom: relative difference 
between MCSM and FCI calculations of the energy (black) 
and rms matter radius (red). 



bling of the MCSM basis dimension leads to a reduction 
of the difference with the FCI results by (approximately) 
a factor of two. However, at A'ghcii — 4 one needs to 
more than double the MCSM basis dimension in order to 
improve the accuracy by a factor of two. 

Naively, one might expect that the number of Monte 
Carlo states needed for a given level of accuracy increases 
with A (and with A'shcii) proportional to the dimension of 
the underlying FCI basis space. However, in practice it 
turns out that the number of required Monte Carlo states 
increases much slower with A than the FCI dimension. 
Note that as A increases, the number of pairwise corre- 
lations grows as and one might expect to require a 
similar increase in Monte Carlo basis states in order to 
achieve a given level of accuracy with strong NN inter- 
actions. Hence one could expect a much more modest 
increase in the number of required Monte Carlo states 
for a given accuracy than the dramatic growth with A of 
the dimension of the underlying FCI basis at fixed iVshoii, 
see Table m Indeed, for A'^ghcii = 3 this dependence seems 
to be roughly between linear and quadratic in A, though 
for A^sheii = 4 the trend is not very clear. Also, so far 



we have only looked at p-shell nuclei, and it is as of yet 
unclear how convergence behaves in the sd shell. 



C. Extrapolation to FCI 

To obtain the converged energy at fixed A^shoU we ex- 
trapolate the MCSM results by using the energy variance, 
which is a new ingredient of the MCSM approach p7| . 
The energy variance AE2 is defined as 



AE2 = {^\H^\^) - {{^\H\^)f 



(6) 



For an eigenstate of H, the energy variance is zero, but if 
^' is not an exact eigenstate of H the energy variance is 
larger than zero. As we increase A^;,, the number of Monte 
Carlo states in the MCSM calculations, we get a better 
and better approximation of the (lowest) eigenstates of 
H. Therefore, AE2 approaches zero from above as A^ 
increases. We use this to obtain an estimate of the exact 
FCI answer. 

We plot the MCSM results for the ground state energy 
of ^Li at different values of A^ as a function of the evalu- 
ated energy variance AE2, see Fig. [71 For AghcU = 2 and 
3 (red and green symbols), the MCSM energy converges 
rapidly to the FCI result (top panel), and the energy 
variance goes to zero (bottom panel). For A^sheii = 4 and 
5 (blue and purple symbols), the energy variance does 
decrease with increasing Nh, but does not reach zero in 
our calculations. For comparison, the open symbols at 
AE2 ~ are the results of our (exact) FCI calculations. 

The behavior of energy as function of the energy vari- 
ance is monotonic and can be extrapolated to zero en- 
ergy variance (which corresponds to the exact e nerg y) 
by quadratic fitting functions as was done in Ref. [17j . 



E{AE2) = Co + C1AE2 + c2{AE2f 



(7) 



with the fit parameters, cq, ci, and C2. Here, cq gives the 
exact energy, E{AE2 = 0). Indeed, the extrapolations 
for NshcW = 4 and 5 reproduce the exact FCI results to 
within a few tens of keV, well within the numerical uncer- 
tainty in the extrapolation. The numerical uncertainty 
for the extrapolation is estimated based on the uncer- 
tainties Sci in each of the three fit parameters Ci of the 
quadratic fit. We treat these three uncertainties as in- 
dependent, and combine them at the MCSM result with 
minimum energy variance, x = min(Ai52), to produce an 
overall estimate of the extrapolation uncertainty 




6cr 



5E_ 

5ci 



5E 

SC2 



(8) 

The FCI and the MCSM results for the energies with 
and without the energy variance extrapolation are all 
summarized in Table IIVI Note that we also quote the 
estimated uncertainty from extrapolation in Table IIVI 

We use a similar extrapolation for the rms matter radii 
and, if possible, also for the magnetic dipole and electric 
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FIG. 5: (Color online) The relative difference between MCSM and FCI calculations of the energy (black), rms matter radius 
(red), and magnetic moment (blue) for several A'shcii values as function of the number of MC basis states, Nb. Top to bottom 
on the left: '^Reir = 0+), '^Li(l+), ^Li(|"), and ''Lid"); Top to bottom on the right: *Be(0+), ^°B(1+), i''B(3+), ^^0(0+). 



quadrupole moments. However, the approach of these 
observables to the exact FCI result is generally not mono- 
tonic, and therefore not as easy to extrapolate. In prac- 
tice we use a linear extrapolation for these observables, 
and apply the extrapolation only if the energy variance 
plot appears to be linear. The detailed dependence of 
both the energy and the other observables on the energy 
variance is presented in the Appendix □ 



D. Point- particle rms radii 



We present the point-nucleon rms matter radii in Fig. [5] 
and Table |V] calculated with the wave functions of the 
MCSM and FCI methods. For this comparison, we eval- 
uate the rms radius of the internal degrees of freedom 
— that is we use the radius operator that depends only 
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FIG. 6: (Color online) Number of MC basis states, A^i,, re- 
quired for a given accuracy of the MCSM energy for (a) 
iVehoii = 3 and (b) iVehoii = 4 for *He(J'^ = 0+), «He(0+), 
«Be(0+), and "-^€{0+). 



on the coordinates with respect to the cm. of the sys- 
tem. Thus, ahhough the nuclear wave functions contain 
mixtures of various components of cm. motion, the use 
of the internal coordinates for the radial operator will 
provide a more accurate rms radius for eventual compar- 
ison with experiment. In addition, at the present level of 
benchmark effort this is sufficient to compare results be- 
tween these approaches. As mentioned above, the exact 
separation of the cm. motion from the internal motion is 
a nontrivial challenge for the MCSM and FCI approaches 
while that separation may be assured in the NCFC ap- 
proach by use of a constraint on the cm. motion. 

The MCSM results in Fig. Hand those in Table Ella- 
belled "extrp" are obtained by extrapolation with first- 
order polynomials using their dependence on the energy 
variance (see the Appendix for more details). We find 
the differences between the extrapolated MCSM and FCI 
rms matter radii to be less than 0.1%, and within the esti- 
mate of the extrapolation uncertainty. As a consequence, 
the open symbols for FCI lie nearly on top of the solid 
symbols for the extrapolated MCSM so that they are not 
separately visible in Fig. HI However, the MCSM results 
for ^He and ^Li in A'shoii = 5 with only 50 Monte Carlo 
basis states are not sufficient for an extrapolation to the 
exact FCI result; more Monte Carlo states are needed 
for a reliable extrapolation for these cases. Note that the 
rms results for -'^^B and ^^C at A^sheii = 4 were obtained 
only within the MCSM approach. 



FIG. 7: (Color online) The convergence of the MCSM ground 
state energy of ^Li to the FCI result; circles (red), triangles 
(green), squares (blue), and diamonds (purple) indicate the 
results at A'shoii = 2, 3, 4 and 5, respectively, (a) Ground 
state energy as function of the number of Monte Carlo basis 
states, A'i,; (b) Energy variance AE2 and extrapolation to 
zero energy variance. 



E. Dipole and quadrupole moments 

In Fig. [9] we plot the MCSM results for the magnetic 
moment of ^Li as function of Ni, (top) and as function of 
the evaluated energy variance Ai?2. For iVshcii = 2 and 
3 (red and green symbols), the MCSM results converge 
rapidly to the FCI result (top panel), and the energy 
variance goes to zero (bottom panel). For A'sheii = 4 
(blue symbols), the MCSM results do seem to converge 
to the FCI results, and with a linear extrapolation on 
the energy variance we get good agreement with the FCI 
results. However, just as for the rms radius, the MCSM 
results for /i at A'shcii = 5 with only 50 Monte Carlo basis 
states are not sufficient for an extrapolation to the exact 
FCI result; more Monte Carlo states are needed for a 
reliable extrapolation for this case. 

We summarize our comparison of MCSM and FCI 
results for the magnetic dipole moments and electric 
quadrupole moments in Table IVll both with and without 
the extrapolations of the MCSM results. We use a lin- 
ear extrapolation using the energy variance, see the Ap- 
pendix for more details. As mentioned above, we cannot 
perform a reliable extrapolation from the MCSM results 
to the exact FCI result for the dipole and quadrupole 
moments of ^Li at A^shoii = 5 with only 50 Monte Carlo 
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FIG. 8: (Color online) Comparisons of the point-particle rms 
matter radii between the MCSM and FCI. Circles (red), tri- 
angles (green) , squares (blue) and diamonds (purple) indicate 
the results at A'^shcii = 2, 3, 4 and 5, respectively; solid (open) 
symbols stand for the MCSM (FCI) resuhs. The MCSM re- 
sults are those after extrapolation as described in the text. 
Note that the '^^B and ^^C results at A^shcii ~ 4 were obtained 
only within the MCSM. 



Nuclei 


Method 


A'^sholl ~ 3 


Afshoii = 4 


A^'shcii = 5 


"He (0+) 


MCSM 


1.355 


1.379 


1.409 




extrp 






1.410(1) 




FCI 


1.355 


1.379 


1.410 


^He (0+) 


MCSM 


1.843 


1.811 


1.864 




extrp 


1.843(1) 


1.813(1) 






FCI 


1.843 


1.813 


1.881 


•^Li (1+) 


MCSM 


1.871 


1.842 


1.889 




extrp 


1.871(1) 


1.846(1) 






FCI 


1.871 


1.846 


1.913 


'Li (1/2-) 


MCSM 


1.958 


1.921 






extrp 


1.959(1) 


1.925(1) 






FCI 


1.959 


1.926 




^Li (3/2-) 


MCSM 


1.931 


1.895 






extrp 


1.932(1) 


1.900(1) 






FCI 


1.932 


1.901 




*Be (0+) 


MCSM 


1.831 


1.958 






extrp 


1.831(1) 


1.960(1) 






FCI 


1.831 


1.960 




i«B (1+) 


MCSM 


1.834 


1.936 






extrp 


1.836(1) 


1.967(2) 






FCI 


1.836 






i«B (3+) 


MCSM 


1.829 


1.909 






extrp 


1.830(1) 


1.926(1) 






FCI 


1.830 






i^C (0+) 


MCSM 


1.722 


1.820 






extrp 


1.723(1) 


1.833(1) 






FCI 


1.723 







TABLE V: Point-particle rms matter radii (in fm) evaluated 
relative to the nuclear cm. The entries labeled "MCSM" in- 
dicate the MCSM results before the energy variance extrapo- 
lation, while the those of the "extrp" line denote the MCSM 
results after the extrapolations. See Table HVl for A*";, and hu 
values. 
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FIG. 9: (Color online) The convergence of the MCSM mag- 
netic moment of ^Li to the FCI result; circles (red), trian- 
gles (green), squares (blue) and diamonds (purple) indicate 
the results at A'^shoii = 2, 3, 4 and 5, respectively, (a) Mag- 
netic moment as function of the number of Monte Carlo basis 
states, Nb; (b) Energy variance and extrapolation to zero en- 
ergy variance. For comparison, the dashed lines (top) and 
open symbols at AE2 = (bottom) are the results of our 
(exact) FCI calculations. 



states. The only other case for which we could not per- 
form a reliable extrapolation is the magnetic moment of 
the (lowest) l"*" state of ^°B. This is likely to be related 
to the fact that there are two 1+ states relatively close to 
each other (experimentally their energies differ by about 
1.5 MeV). 

We again find the differences between the MCSM and 
the exact FCI results to be small, typically 1% or less 
with 100 Monte Carlo states, both for the magnetic 
dipole moments and for the electric quadrupole mo- 
ments. The exception is the quadrupole moment of the 
ground state of ^Li: at A'shcii = 5, the difference be- 
tween the MCSM and FCI calculations is almost 10% 
with 50 Monte Carlo basis states. Note however that this 
quadrupole moment is exceptionally small in magnitude: 
although the relative difference between the NCSM and 
the FCI result is significantly larger than for most other 
observables we have looked at, the absolute difference is 
rather small. 

A linear extrapolation using the energy variance brings 
the MCSM results even closer to the FCI results, see the 
Appendix for more details. However, the results after 50 
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Nuclei 


Method 


-^shcll — 3 


A^shell = 4 


A''shcii = 5 


"Li (1+) 


MCSM 


0.836 


0.834 


0.836 




extrp 


0.835(1) 


0.833(1) 






FCI 


0.835 


0.832 


0.830 


•^Li (1/2-) 


MCSM 


-0.842 


-0.816 






extrp 


-0.840(1) 


-0.806(2) 






FCI 


-0.840 


-0.807 




^Li (3/2-) 


MCSM 


3.061 


3.025 






extrp 


3.057(1) 


2.995(2) 






FCI 


3.056 


2.993 




(1+) 


MCSM 


0.503 


0.533 






extrp 


0.508(1) 








FCI 


0.509 






i«B (3+) 


MCSM 


1.820 


1.814 






extrp 


1.818(1) 


1.819(1) 






FCI 


1.818 






Q (efm^) 


"Li (1-^) 


MCSM 


-0.259 


-0.282 


-0.276 




extrp 


-0.259(1) 


-0.285(1) 






FCI 


-0.259 


-0.285 


-0.302 


^Li (3/2-) 


MCSM 


-1.766 


-2.006 






extrp 


-1.750(1) 


-1.958(3) 






FCI 


-1.750 


-1.940 




i«B (1+) 


MCSM 


-1.712 


-2.417 






extrp 


-1.703(2) 


-2.436(8) 






FCI 


-1.698 






(3+) 


MCSM 


3.532 


5.222 






extrp 


3.503(1) 


5.250(11) 






FCI 


3.503 







TABLE VI: Dipole (top) and quadrupole moments (bot- 
tom) calculated using the wave functions obtained within the 
MCSM and FCI methods. The entries of the MCSM indicate 
the MCSM results before the energy variance extrapolation, 
while the those of the "extrp" line denote the MCSM results 
after the extrapolations. See Table [TVl for Ni, and hw values. 

Monte Carlo states for ^Li at A^shoii = 5 are not suffi- 
ciently close to the FCI result to do such an extrapola- 
tion. 

The magnetic dipole moments, top panel of Fig. I10| 
depend only very weakly on the basis space trunca- 
tion parameters A'shcii and hw — much weaker than the 
quadrupole moments, bottom panel of Fig. IIOI and than 
the rms radii. Fig. [8l The differences are less than 2% 
and are not visible on this scale. Furthermore, the dipole 
moments are in very good agreement with the NCFC re- 
sults, which means that they are converged to within a 
few percent with respect to the basis space truncation. 
Our results with JISP16 are also in good agreement with 
the available data for the ground states. 

On the other hand the quadrupole moments do de- 
pend significantly on the basis space truncation parame- 
ters A'shcii and hio, as can be seen from the bottom panel 
of Fig.lTUland from Table IVTl This is not surprising, given 
the dependence of the rms radius on the truncation pa- 
rameters, and given the fact that the quadrupole moment 
receives significant contributions from the asymptotic tail 
of the wave function, which is not very well represented 
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FIG. 10: (Color online) Comparisons between the extrap- 
olated MCSM and FCI results for (a) the magnetic dipole 
moments, and (b) electric quadrupole moments. The conven- 
tions for the symbols are same as in Fig. [S] crosses indicate 
the experimental values for the ground states from Ref. [2^ . 



in a HO basis. One needs to include (much) higher HO 
shells in order to build up a realistic tail for the wave 
functions. Nevertheless, our results are in qualitative 
agreement with the available data: small and negative 
for the ground state of ^Li, large and negative for the 
ground state of ^Li, large and positive for the ground 
state of i"B. 



V. SUMMARY AND OUTLOOK 

We have performed benchmark calculations of the en- 
ergies, point-particle rms matter radii, and electromag- 
netic moments for nine states in light nuclei ranging from 
''He to ^^C. Where possible, we have solved for these 
properties using the FCI, MCSM and NCFC approaches. 
The energies and the point-particle rms matter radii cal- 
culated by MCSM were extrapolated as a function of 
energy variance. All results are found to be consistent 
with each other to within quoted uncertainties when they 
could be quantified. Where we could not obtain quanti- 
fied uncertainties, the results were found to differ typi- 
cally by a few percent among the available methods with 
very few exceptions. The MCSM and FCI results are 
very close to each other with small differences (of a few 
percent in most cases) arising mainly from the limited 
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number of MCSM basis sampled stochastically for diag- 
onalization and from MCSM energy variance extrapola- 
tion uncertainties. We include selected NCFC results in 
order to gauge the increases in basis spaces needed to 
better approach the fully converged results (basis space 
cutoff independence) in future efforts. 

Since the MCSM computational effort scales more fa- 
vorably with increasing basis space and increasing nu- 
cleon number, we expect that the MCSM will further 
develop into a powerful tool for ab initio nuclear theory. 
To reach this goal, we will need to expand the basis space, 
treat the role of cm. motion and include the Coulomb in- 
teraction as well as NNN interactions. These challenges 
will be addressed in future efforts. 
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Appendix: Extrapolation of MCSM results 

With increasing Monte Carlo basis dimension Nb, the 
MCSM results converge to the FCI results. In order 
to obtain an estimate of that exact FCI answer, we ex- 
trapolate the energy and other observables evaluated by 
MCSM using the energy variance. That is, the MCSM 
results are plotted as a function of the evaluated energy 
variance 

A£;2 = - ((*|iJ|*))^ (A.l) 

and then extrapolated to zero variance as we show below. 
We also investigate the uncertainties of this extrapolation 
and report those uncertainties in Tables IIVHVII 



Figure [TT] shows the energies as function of the energy 
variance for ^He, ^Li, "^Li, ^Be, ^°B, and ^^C. For iVshoii = 
2 and 3 there is no need for any extrapolation: with 
100 Monte Carlo states, there is very good agreement 
between the MCSM results and the FCI results. For 
A'shoU = 4 and 5 we use a quadratic polynomial fit to 
extrapolate AE2 to zero. We also make an estimate of 
the numerical uncertainty in this extrapolation. These 
extrapolated MCSM results are in good agreement with 
the available FCI results (indicated by the open symbols 
at AE2 =0). In Table |IV] we give both the MCSM 
results, and the extrapolated results with extrapolation 
uncertainty. 

Figure [T^ shows the rms matter radii as function of the 
energy variance for ^He, ^Li, ^Li, *Be, ^"B, and ^^C. For 
A'shcii = 2 and 3 there is no need for any extrapolation: 
with 100 Monte Carlo states, there is very good agree- 
ment between the MCSM results and the FCI results. For 
A'shcU = 4 we use a linear fit to extrapolate AE2 to zero. 
We also make an estimate of the numerical uncertainty 
in this extrapolation. These extrapolated MCSM results 
are in good agreement with the available FCI results (in- 
dicated by the open symbols at AE2 =0). In Table IVl 
we give both the MCSM results, and the extrapolated 
results with extrapolation uncertainty. 

Unfortunately, 50 Monte Carlo states is not sufficient 
to extrapolate the radii of ^Li and ^He at A^sheii = 5: 
the purple diamonds in the upper left figures cannot be 
extrapolated reliably to AE2 = 0. This is also the case 
for the magnetic dipole moment, see Fig. 1101 and the 
electric quadrupole moment, see Fig. 1131 below. 

Figure [13] shows the electric quadrupole moments as 
function of the energy variance for the states that have 
J > 1. For A'shcii — 2 there is no need for any extrapo- 
lation. However, both the A'^jhoii — 3 and the A^'sheii = 4 
MCSM results with 100 Monte Carlo states can be im- 
proved by a linear extrapolation. As already mentioned, 
50 Monte Carlo states is not sufficient to extrapolate the 
quadrupole moment of ^Li at A^shcU — 5. 

Finally, in Fig. [14] we show the magnetic dipole mo- 
ments as function of the energy variance for the ^Li and 
-^°B states that have J > 1. Again, for A^shcii = 2 there is 
no need for any extrapolation. Both the A^shoii — 3 and 
the TVshcii = 4 MCSM results with 100 Monte Carlo states 
can be improved by a linear extrapolation. However, the 
dependence of the magnetic moment of the (lowest) 1+ 
state of ^"B does not seem to converge as the energy vari- 
ance decreases. This is possibly caused by the proximity 
of a second 1+. 

The extrapolated MCSM results for both the magnetic 
moments and the quadrupole moments are in good agree- 
ment with the available FCI results (indicated by the 
open symbols at AE2 — 0)- In Table IVTl we give both 
the MCSM results, and the extrapolated results with ex- 
trapolation uncertainty. 
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FIG. 11: (Color online) The energy variance and extrapolation to the FCI result for the energies. Circles (red), triangles 
(green), squares (blue) and diamonds (purple) indicate the results at A'shcii = 2, 3, 4 and 5, respectively; the open symbols at 
AE2 = are the exact FCI energies. Top to bottom on the left: ^lle{J'^ = 0+), ®Li(l+), '^'L\{\~), and ''Li(|"); Top to bottom 
on the right: *Be(0+), '°B(1+), i''B(3+), and ^^C{0+). 
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FIG. 12: (Color online) The energy variance and extrapolation to the FCI result for the rms radii. Symbols are the same as 
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FIG. 13: (Color online) The energy variance and extrapolation to the FCI result for the quadrupole moments. Symbols are the 
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;i") and ^Li(f 



as in Fig.[n] Top and bottom on the left: ■^Li(i ) and ^Li(| ); Top and bottom on the right: ^''B(l+) and ^°B(3+) 



19 



[1] S. C. Pieper, R. B. Wiringa, and J. Carlson, Phys. 
Rev. C 70, 054325 (2004); K. M. NoUett, S. C. Pieper, 
R. B. Wiringa, J. Carlson, and G. M. Hale Phys. Rev. 
Lett. 99, 022502 (2007); S. C. Pieper, in Proceedings 
of the International School of Physics "Enrico Fermi", 
Course CLXIX, edited by A. Covello, F. lachello and 
R. A. Ricci (Societ Italiana di Fisica, Bologna, 2008), 
p. Ill, arXiv:071 1.1500 [nucl-th]; reprinted in La Riv- 
ista del Nuovo Cimento, 31, 709, (2008), and references 
therein. 

[2] P. Navratil, J. P. Vary, and B. R. Barrett, 
Phys. Rev. Lett. 84, 5728 (2000); Phys. Rev. C 
62, 054311 (2000); S. Quaglioni and P. Navratil, 
Phys. Rev. Lett. 101, 092501 (2008); Phys. Rev. C 79, 
044606 (2009). 

[3] G. Hagen, T. Papenbrock and M. Hjorth- Jensen, Phys. 
Rev. Lett. 104 182501(2010), and references therein. 

[4] E. Epelbaum, W. Glockle, and Ulf-G. Meissner, Nucl. 
Phys. A 637, 107 (1998); 671, 295 (2000). 

[5] D. R. Entem and R. Machleidt, Phys. Rev. C 68, 
041001(R) (2003). 

[6] R. B. Wiringa, V. G. J. Stoks and R. Schiavilla, Phys. 
Rev. C 51, 38 (1995). 

[7] S. C. Pieper, V. R. Pandharipande, R. B. Wiringa, and 
J. Carlson, Phys. Rev. C 64, 014001 (2001) 

[8] S. C. Pieper, AIP Conf. Proc. No. 1011, 143(2008). 

[9] A. M. Shirokov, J. P. Vary, A. I. Mazur and T. A. We- 
ber, Phys. Letts. B 644, 33 (2007); A. M. Shirokov, 
J. P. Vary, A. I. Mazur, S. A. Zaytsev and T. A. We- 
ber, ibid. 621, 96 (2005); subroutines to generate this 
interaction in the relative-center-of-mass HO basis are 
available at nuclear.physics.iastate.edu. 
[10] M. Honma, T. Mizusaki and T. Otsuka, Phys. Rev. 
Lett. 75, 1284 (1995); 77, 3315 (1996); T. Otsuka, 
M. Honma and T. Mizusaki, ibid. 81, 1588 (1998); 
for review and further references, see T. Otsuka, 
M. Honma, T. Mizusaki, N. Shimizu, and Y. Utsuno, 
Prog. Part. Nucl. Phys. 47, 319 (2001). 



[11] P. Maris, J. P. Vary, A. M. Shirokov, Phys. Rev. C 79, 

014308 (2009); P. Maris, A. M. Shirokov and J. P. Vary, 

ibid. 81, 021301(R) (2010); C. Cockrell, J. P. Vary and 

P. Maris, ibid. 86, 034325 (2012). 
[12] G. Hagen, T. Papenbrock and D. J. Dean, 

Phys. Rev. Lett. 103, 062503 (2009). 
[13] J. P. Vary, The Many Fermion Dynamics Shell Model 

Code, Iowa State University, 1992 (unpublished); 

J. P. Vary and D. C. Zheng, ibid., 1994 (unpublished); 

P. Sternberg, E. G. Ng, C. Yang, P. Maris, J. P. Vary, 

M. Sosonkina, and H. V. Le, in Proceedings of the 2008 

ACM/IEEE conference on Supercomputing (IEEE Press, 

Piscataway, NJ, 2008), pp. 15:1-15:12. 
[14] P. Maris, M. Sosonkina, J. P. Vary, E. G. Ng and C. Yang, 

International Conference on Computer Science, ICCS 

2010, Procedia Computer Science 1, 97 (2010). 
[15] N. Shimizu, Y. Utsuno, T. Abe, and T. Otsuka, RIKEN 

Accel. Prog. Rep. 43, 46 (2010). 
[16] Y. Utsuno, N. Shimizu, T. Otsuka and T. Abe, Com- 

put. Phys. Comm. 184, 102 (2013). 
[17] N. Shimizu, Y. Utsuno, T. Mizusaki, T. Otsuka, T. Abe, 

and M. Honma, Phys. Rev. C82, 061305(R) (2010); AIP 

Conf. Proc. No. 1355, 138 (2011). 
[18] T. Abe, P. Maris, T. Otsuka, N. Shimizu, Y. Utsuno, and 

J. P. Vary, AIP Conf . Proc. 1355, 173 (2011). 
[19] G. Puddu. [aTXr^TT201 .0600 [nucl-th]. 
[20] R. Roth, Phys. Rev. C 79, 064324 (2009). 
[21] T. Dytrych, K. D. Sviratcheva, C. Bahri, J. P. Draayer 

and J. P. Vary, Phys. Rev. Lett. 98, 162503 (2007); 

T. Dytrych, K. D. Sviratcheva, C. Bahri, J. P. Draayer 

and J. P. Vary, J. Phys. G. 35, 095101 (2008); ibid. 35, 

123101 (2008). 

[22] L. Liu, T. Otsuka, N. Shimizu, Y. Utsuno and R. Roth, 

Phys. Rev. C 86, 014302 (2012). 
[23] N. J. Stone, Atomic Data and Nuclear Data Tables 90, 

75 (2005). 



